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Locatjon  of  sales  outlets  Is  of  major  concnrn  to 
such  orp,  rTnl>!ations  as  oil  companies,  banks,  etc.  The  problem  Ir,  of- 
ten approached  in  two  steps.  The  total  r.^irkct  is  divided  into  r<.>- 
p, ions,  and  the  first  stape  amounts  to  decidinp  in  which  repions  to 

expand  or  to  contract.  Thus  for  example  in  period  t,  the  company 

1 1  2  t 

(i)  '.•;ill  add  n.     outlets  in  area  1,  n^   in  area  2,  etc.  We  will 
1  i 

refer  to  the  first  step  as  the  ap  pre  pate  location  prohlern.  The  se- 
cond stape  consists  in  chooLlnp  specific  sites  for  these  new  cut- 
lets. In  practice,  these  decisions  are  rep,  arded  as  rather  in  da  pen- 
dent, especially  because  they  are  macle  at  different  levels  in  the 
orpani  ;^a  t  i  on  .  The  number  of  new  outlets  is  a  corporate  decision, 
whereas  specific  sites  are  selected  at  the  repional  level,  subject 
however  to  approval  by  the  corporate  headnucrtsrs.  Many  companies 
feel  that,  at  least  for  the  time  beinp,  hiierarchical  linFing  of  the 

aggregate  and  detailed  problerris  is  neither  worth  the  effort  ncr  the 

^  1 

cost  . 

In  this  paper  our  sole  crncern  will  be  with  aggi'cgs- 
te  location.  Our  procedure  will  be  closely  related  t'j  a  model  deve- 
loped by  Hartung  and  Fisher  f/'].  Hov-jcver,  their  work  lacked  robust- 
ness, and  suffered  from  a  variety  of  deficiencies  :n  the  estimation 
of  the  model  parameters.  In  section  2    ^.'e  will  review  the  Hartunp- 
Fisher  [hereafter  M-F]  model  and  the  various  weaknesses  associated 
v^/ith  it.  In  section  3  the  estimation  problems  wjll  be  examined.  In 
section  1  vie    will  propose  various  changes  to  the  r'^dsl  whicti  v.' ill 
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maKc  it  robust,  ond  we  will  use  data  from  a  major  oil  corpany  in 
a  European  country  to  est  i  ma  to  the  parr.mGtcrr.  and  validate  cur  ap- 
proach . 

2 

^  •  THE  HARTDNH-FISHrR  rODFL 

The  market  is  reductrd  to  a  quasi  duopoly,  that  is, 
we  consider  our  brand.l,  versus  competitive  brands  taken  top, ether, 
c.  Buyinr  behavior  is  described  as  a  first  order  Markov  chain. 
The  transition  prnbabilit,  ins  are  defined  as  follows  : 

X.  =  probability  that  a  person  i-;ho  usual  I  y     buys  brand  i 

in  period  t-1,  will  usual  7- u    buy  brand  i  in  period  t. 

Adding  the  word  "usually"  broadens  the  definition  used  by  H-F  be- 
cause it  alio  v;  3  for  incidental  purchases  of  a  competitive  brand. 
-This  is  important  for  a  product  such  as  gasoline.  Take  a  person 
v;ho  usually  buys  brand  1.  On  a  given  day  he  is  running  cLit  of  gas 
on  a  thruway  and  tanks  at  the  next  s'ervice  area.  If  the  brand  is 
not  i,  we  should  not  conclude  yet  that  brand  switching  has  occurred 


The  other  transition  probabilities  are  similarly 


defined  : 


probability  of  switching  from  c  to  i 

probability  of  remaining  a  buyer  of  a  competitive  bran r I 

pro  [j  ability  of  switching  from  i  to  c. 


Market  share  of  bronri  1  in  period  t  Is 


1)     m=x.*m_+o.-m 

i.t     1    i,t-1      1    c,t-1 


In  steady  state,   rn 


i ,t  i  .t-1 


,     and     thus 


(2)  m  =    a  .     /     (1     -     X  .     +    o  .  )     =    c^.     /     fo       ♦    a     ) 

i,ei  11  1  ic 


The     transition    probabilities     v-/ill     be     functions     of 
the    decision    variables    of    company    i.     such    as    advertising,    cxpendi- 


.examp 1 e , 


.     =    s.     (a.,     d.,     p.,     a     ,     d     ,     p 
i  1  i  I'^ic  cc 


H-F  consider  only  the  number  of  sales  outlets  as  determinants  of 
the  transition  probabilities.  The  following  functions  v;ere  post  u  la- 


(3)   ^i  =  1^1  •  d^  /  (d^  *  ^i^ 


d  .  /  [d   +  d  .  ) 

1        c       1 


5. 


For  a  given  value  of  d  ,  both  X   and  a.  are  increasing  with  ri, 

and  show  decreasing  returns.  The  limits  of  X.  and  of  o.  for  d   ->•  <» 

i  11 

are  k   and  K   respectively.  However,  K   and  K  ,  and  hence  X.  and 
a.     are  not  restricted  to  he  betv^een  zero  and  one.  For  example,  H-F 


larger  than  one.  Since  X   is  a  probability,  its  dependence  on  d. 

and  d   should  be  constrained  in  such  a  way  that  its  value  will  lie 
c 

in  the  [0,1]  interval.  The  H-F  model  therefore  loc^s  robustness. 
H-F  are  well  aware  of  that.  They  state  that, 

"Unless  k   =  k   =  1.0,  X.  and  o,  are  not  probabilities  for  all 

values  of  d./(d   +  d,].  Ho v; ever,  equations  '3)  and  (4}  can  be 
1    c     i 


ted. 


More  general  functions  can  be  substituted  for  equations  (3)  and 

(4)  without  invalidating  later  results.  The  authors  found  that 

3 

for  their  problem  equations  (4]  were  sufficiently  accurate". 


Section  4  v-jill  be  devoted  to  examining  ways  through  v-.  ch   the  model 
can  be  mads  robust. 

Let  us  now  relate  market  share  to  the  parameters 


k   and  k.  Substituting  {3}  and  (4]  for  X   and  o   in  ['^)     gives 


6. 


(5) 


i.e 


C1 


Lpt  q  -    sales  of  i,  and  Q  =  industry  sales.  Market  share   m.  =  q./ 

Note  that  H-F  replac::'   m^     in  (51  by  a./O.  We  should  observr;  that 
I.e  -^1 

this  implies  the  assumption  that  observed  market  share  is  equal  to 
steady  state  market  share.  V-'e  will  return  to  this  issue  in  ssctiori  • 


terms,  H-F  obtain. 


(6] 


q^/d^  =  k^  •  Q  /  (d^ 


From    (65     it    follows    that    for    a    piven    value    of    d     ,     and    v/ithi    p    <    D, 

'-  c 

average  sal  as  per  outlet  will  increase  with  ri,  and  f.  oes  to  infini- 
ty when  (d^  +  y«d.)  goes  to  ^ero.  With  k   =  4.44  and  k   =  0.64, 
p  =  -2.a0,  anc)  therefore,  the  model>woulri  predict  infinite  sales 

per  outlet  when  d./d   =  1/2.8.  For  values  d  . /d   >  1/2.8,  the  mo- 
^  1   c  1   c 

del  would  predict  negative  sales.  The  function  pattern  is  depicted 
inFigurol. 


Insert  Fioure  1  cbD'jt  here 

On    the    other    hand    i f    p    >    0 ,     the    model     average    sal^' 
pf  1     -^'itlet    will     decrease    when    tlic    nufitaer    of    outlets     incrftastrG. 
l.ambin    applied     u'-"    H-F    model    on    a    brand    of    p;asoline    in    a    European 


7. 


cruntrv  [10,  chanter  7^,     and  founci  values  of  k 


and  y  =  0.9  5  7.  However,  even  with  p   >  0,  one  may  still  run  into 
difficulties.  With  d.  larp,c  (for  a  given  d  },  q.  approaches  k*  Q/u 


sales  exceed  industry  sales. 

The  next  issue  deals  with  the  estimation  of  the  pa- 
rameters k   and  k„.  Equation  (B)  is  nonlinear  In  K   and  k,  and 
could  not  be  estimated  by  linear  regression.  H-F  applied  a  series 
of   trans  for  rr.ations  to  [6)  and  ended  up  with  the  folio  i-j  in  rr  function 

(7) 


H-F  use  linear  regression  to  estimate  a  and  B.  The  paramete-s  k 
and  k_  are  then  uniquely  determined  from  a  and  B.  There  arf  various 
probleins  absociated  witii  the  Gj.LimdLion  pioceduie  usyd  ij y  H-F. 
These  will  be  examined  in  section  3. 


in  an  aggregate  retail  outlet  location  model  in  i^jhich  the  ob, 
1^  to  maximize  discounted  return  of  firm  (i). 


Let   J     =  the  number  of  regions  (J=1,...,J1 

T     =  the  tine  horizon       (t-0....,T] 

T     =  the  disco  LI  nt  factor 

r;:    =  return  per  unit  sold  in  rep, ion  J,  period  t 

ri'T  ,  q.  ,  ...  as  d.,  q.  but  superscripted  for  regions  nnd  time 

d.    =  existing  r.  upher  of  outlets  at  time  D 


Jt 


Jt 


.Jt 


and  in  period  t 


The  return  maximization  model  is  then 


Max  Z    I 
J  t 


.f    ^i'   .f 


subject  to 


■it   jt 
I       c~        n^   <  b . 

J    '    ^   "   ' 


qf  =  K^  .  qJ^  /  [dJ^  .  ..dfl 


Vj.t 


df  =  dJ°  .   I       nJ^    Vj.t 


Jt  ^   jt 
'i   -   i 


Vj.t 


9. 


The  readRr  wil]  notice  that  with  p  <  D,  n-j       is  an 

'•it  it 

increasing  function   of  d  .   and  therefore  unle^js  r;:   is  carefully 

specified  as  a  net     return  and  unless  the  discounted  construction 

costs,  i.e.  L    Z     \- c.   n^  ,  are  deduced,  the  objective  func- 

j  t  i^  *  ^J    1    i. 

tion  is  monotcjne  increasing  and  the  solution  is  rrerely  p,  i\'en  by 
the  constraints.  It  is  also  important  to  note  th^t  even  when  the 
discounted  coiistruction  costs  ore  deduced,   the  objective  function 
may  still  be  monotone  increoslnp, ,  for  d^   <   |d  /  \i\  • 


3.  ISSUES  IN  ESTIMATING  THE  H-F  MfinEL 


timated  values  of  a  and  B  in  equation  (7) 


rif;;ht  hand  side  of  the  equation  while  in  fact  the  number  of  retail 
outletsis  expected  to  have  a  causal  effect  on  market  share  and  not 


left  hand  side  of  equation  (7)  is  q./'J  =  (q,/Q)(1  +  d  /d.}.  Or.  equ 

1         1  c   1 


tion  (7)  can  be  rewritten  as 


From  a  causality  point  of  viev^/  we  would  like  to 

have  ri./d   on  the  ripht  and  m.  on  the  I'^ft  of  the  equation.  Scne 
1   c  "^  1 

simple  manipulations  of  equation  (6)  result  in  the  followinp,  lineai 
equation, 


10, 


(9)       1/m.=Y*<5'(cl/ri.) 
1  c   1 


where  y    =     (1  -  P)/nc,  and  6  =  1/ci.  Or  in  ternG  of  k   and  K,  K=  1/ 


the  market  share  function,  wt.'  can  estimate  it  directly  as  nbtainer; 
In  equation  (5)  (i-/:th   m,  r-eplacinp,   m.    }  hv  nonlinear  er.  tinatjnn 


ned  fron  equations  (7),  (9],  and  (5K  The  estimations  were  pcrforme 

5 

on  the  TROLL  systpn  ,  The  estimation  proceriurG  in  TROLL  is  br-seri 

on  r^arquardt's  alporithni  for  least-squares  estimation  of  nonlinear 
paramet  ers  [12  1. 

Wg  estimated  eauations  (7),  (91,  and  (5)  for  a  major 
brand  of  gasoline  in  a  Furopcan  coun'try-  There  arc  35  qu^Tterly 
observations,  from  the  first  quarter  of  1962  to  the  third  quarter 
of  1970.  The  statistical  results  for  the  nonlinear  model  are  based 

on  its  linearized  form  around  the  ontimu'n,  i.e.  the  minimum  of  the 

6 

residual  sum  of  squares  .  Let   0  =  (0,  ...  C   ...  0„)   ^  vnctcr  of 


Pt 
The  general  model  is  written  as. 


Let     e     bG     the    final     InrDst-squnr'RS     er-  tirpatG    of     Q  ,     Vlith     0     close     to     3 


9f (y,  ,ei 

(in)        E(\'    }    =    f(y^,o]    ♦    T    \—-~-A -    ^^      -    0 

p    1  p      I    e=G       ^         ^ 


with  f 


pt 


3f (X^.O) 


f^      =   f  (x^  .0)    -    j:     f    ,      0 
t  t  pt        p 

p 


jqufltion     (10)     can    be    v.Titten    as    a    linear    function    of    the    parametBrs 


t  P  pt  t 


wherB,  z   f  V   -  f.   •  The  statistics  for  the  nonlinear  model  are 
similar  to  those  of  a  lineai-  rep;ross'ion.  Thus,  if  we  definr;  the 
following  (T  X  P)  -  matrix  : 


P1 


12     ?? 


/here  T  is  the  number  of  observations  available,  i  f-,  -  n  the  e  s  1 3  r  a  t  e  r 

2-1  7 

^ariance-covariancc    n>atrix     of     0    is     V  [  0)     =     s      •    (F'F)      ,     in    v-zhich     : 


is    tbe    residual    mean    sqLiare 


1?. 


Note  that  n   ir,  no  Icn^Gr   an  unhiased  c.i  timet  n  of 
a^  ,     thR  disturbance  variance  and  that  even  when  the  en- or  tf;rm  f: 
is  normally  diGtributnd,  0  is  no  Innp,  er  normally  dlntrihut.  ed.  As  a 
result,  the  u'.  nal  t-,  F-,  and  Durbin -Watson-  statistics  are  not 
valid  in  psner.-^l.  Howpver,  these  statistics  will  be  reported  here 
they  shoLilri  therefore  he  rep,  arded  as  mere  com.  pa  risen  valur^s. 

The  rf;  suits  for  tirand  i  are  presented  in  Taiile  I. 


share  of  brand  :  to  less  than  5  per  cent.  Yet,  ovpr    the  i%'hole  pe- 
riod of  observations,  its  outlet  share  was  hip  her  than  8  per  cent. 
Worse  even,  k   is  nep,  otive.  The  values  of  the  '"'urbin-VJatson  sta- 
tistic indicates  autncon- elation  of  the  residuals,  and  hence  pro- 
blems with  the  model  specification.  There  are  several  possifile 
explanations.  One  is  that  im,  port  ant  additional  explanatory  varia- 
bles may  have  been  left  out.  A  previous  study  of  these  data  shows 

8 
that  this  is  not  the  case  .  More  liRely  reasons  are  the  misspecifir 

tion  of  the  transition  probability  functions,  and  the  assumption 

that  observed  market  shares  are  enuilibrium  values. 


In  other  cases  one  may  have  better  luck.  For  exam- 

9 
pie,  VIP.    applied  the  H-F  model  on  another  brand  (a)  ,  and  found  re- 
sults similar  to  those  in  the  H-F  paper  as  illustrated  in  Table  II 

The  Durbin- V/atson  statistic,  ho v; over,  ap, ain  indicates  significant 

10 
autocorrelation 


Insert  Tables    I    ?.nd    il    aiiojt    here 

In  the  next  s  c  c  t  i  o  n  '  w  e  w  j  1  1  r  x  p  1  n  r  b  t  u'  o  w  ri  v  s  in 
which  thR  model  can  \r  c    made  robust.  Wo  will  apply  both  these  rridrl 
formulations  to  the  data  for  brand  i. 

4  ,  ALT  E;_P  h!_AT  TVF  FORniJt  -VrrO^ N  S 

In  this  section  we  will  redefine  the  transition 
prohahility  functions.  First,  as   oxponential  functions   of  the 
relative  number  nf  outlets  in  section  4.1.,  and  next  as    lopistic 
f  u  n  c  t  i  o  n  s  j  n  section  4.2.   How  to  rn  ake  use  of  these  models  will  he 
further  examined  in  section  4,3. 

4.1.  Exponential  model 


Define  A .  as 


i 
d./d   approaches  infinity,  A.  approaches  one.  Fnuation  (1?)  thus 

relates  the  relative  number  of  outlets  D.,  to  the  transition  protia- 

bility  X . ,     in  a  robust  way. 

Similarly,  X   is  defined  as 


(13}    \        =  1  -  EXP(-a  •  D  ] 
c  c    c 


0.  And  therefore.   o   =  EXPC-a.  'D.)  ,  a  _,     =    EXPf-a 
c  1    1       i  c 


14  . 


TarKet    share    at    tirrie    t    ic    related    to    market    share 


at    timp    t-1. 


(14)  rn .      .      =     [1     -     E  XP  (  -  a  .    •   n  .     ,  ]     -     E  XP  [  -  a      •   P       ,  )  1    •    m  .     ,     , 

i.t  1  i.t  c  c.t  i,t-1 


+     EXP(- J       •   D       ^ 
c  c.t 


Assume  now  for  a  moment,  as  H-F  did,  that  ohsorvcd 
markot  share  values  pre  Gquilibriurt  values,  for  piven  values  oT 


D.  .  and  D  ^  .     That  is, 
I.t        c.t 


(15) 


1.  ^    =  EXP{-a   •  D   ^]  /  fEXP(-a,  •  D.  ,)  +  EXP(-a   •  D   ^}1 
i.t.e  c    c.t  i    I.t  c    c.t 


Applying  a  logjt  transformation  to  (15)  results  in 


a  linear  mod  e 1   . 


(1G)      l0R[m.  ^    /  (1  -  m.   .    )]  =  a  ,•  D.  ,  -  a   •  n   . 
'^'x.t.G  i.t.e-*     i    I.t     c    c.t 


The  results  of  the  estimation  of 

ted  in  Table  III.  The  estimated  a.  is  negative  which  would  seum    to 

indicate  that  the  assumption  on  enuilibrium  is  not  at  all  Wr-irrcjntec 


It  should  be  c]Gar  that  many  assumptions  -  not 
satisfied  in  this  case  -  have  to  be  made  in  order  to  accept  tf-  e 
equilitjriurn  form  (16).  For  ex  a  n^  pie,  provided  that 

a)  -  the  ccrsumption  patterns  arc  adequately  stable. 


b)  -  the  unit-tirnj  period  is  ouffiniently  long  so  that  the  disrup- 
tive effect -J  of  a  marketing  campaign    launched  in  period 


[t-1,t]  -  on  the  B  teady- st  nt  <;  narknt  shares  can  resorh  In 
a  new  equilibrJum  achieved  within  thq  same  period  ft-l.tl, 

c)     -     the  unit-tirnc  period  is  short  ennuph  so  that  no  compntitive 
reaction  can  interfere  witfi  this  new  ccuilihrium, 

we  may  retain  equation  (15).  The  restrictive  and  somewhat  contra- 
dictory chnr-acter  of  this  no  n -exhaustive  ret  of  assumptions  cxplciin 
why  vm    should  turn  to  dynamic  forms.  This  does  not  imply  that  we 
have  to  disrecnrd  the  steady-state  aspects  w/hen  we  are  about  to 
make  decisions  on  where  to  add  new  outlets.  The  results  in  Tahle 
III  merely  indicate  the  market  dynamics  should  be  taken  into  ac- 
count in  estimatinp,  the  parameters. 

Insert  Table  I  I  I  about  her3 


Equation  (14)  is  intrinsically   nonlinear,  and  was 
.estimated  in  two  different  ways.  First,  we  used  TF^DLL.  Secondly, 

we  applied  the  Senuentia]  Unconstrained  MinJmizotinn  Technique 

1  2 

(SUMT)   .  rjonlincar  programming  can  bn    used  for  nonlinear  estima- 
tion in  the  following  way.  Let  the  model  be. 


r-i  i  n  i  m  i  z  : 


is  achieved  by  solving  the  nonlinear  programming  problem  belcw 


(17)       min   i:J^![   cj(e) 


t  =  1  , 


;rt  Table  IV  and  Fi 


Tnhl'.'  TV  stiov-.'S  twe    re5u]t;;  of  the  n;;  1 1  ma  t  i  on  of 

>quation  (Ih]  usinj;  TROLL  and  SU^iT.  Figure  2  illu?.  tratcs  \.     and  X 

1        c 


VJith  tht:  currRnt  number  of  outlets,  D^ 
A  .  = 


ThGSG  valuR'5  are  very  hip, h,  which  is  to  be  Rxpectcd 

dieted  equilibrium  market  share   is     .  A  one  per  cent  increase 

in  d./Cdj  +  ri  ],  and  assuminr  ri   remains  constant,  equilibrium  mar- 
1    i     c  '   c 

Ket  share  would  increase  by      per  cent.  Whether  such  an  increase 
in  share  would  be  worth v-jile  depends  on  industry  sales  volume,  unit 
profit,  and  the  number  of  new  outlets  needed  to  incr-easc  outlet 
share  tjy  one  percent. 

More  interesting  of  course  is  to  look  at  the  problem 
on  a  regional  basis.  For  this  particular  product,  we  have  informa- 
tion on  four  different  regions.  The  company's  outlet  share  varies 
from  a  high  of  about  ID  per  cent  in  one  region  to  a  low  of  about 
5  per  cent  in  another  region.  Incremental  year'ly  regional  sales 
per  outlet  added  are  sho'.'vn  in  Table  V  for  each  of  the  four  regions. 
Adding  outlets  in  region     contributes  thehighest  marginal  return. 


It  is  quite  possible  that  the  response  parameters 
differ  r  ^  c r o s  s  regions.  For  example,  r u  r- a  1  areas  might  be  d  i  :;  t  i  n - 
gulshed  from  met  rope ] i t an  areas.  In  our  particular  application  in- 


17. 


formation  by  region  wos  .. vailable  only  on  an  annual  basis  over  a 
period  of  six  years.  This  was  insufficient  for  the  purpose  of  es- 
timating responsR  coefficients  by  re  f,  ion. 


Insert  Table  V  about  h^re 


4.2.  L^O£ist_i  c__modRj 

Market  share  as  a  function  of  relative  number  of 

outlets  is  often  thouc^t  of  as  having  an  S  shaped  form.  Various 

oil  companies,  for  example,  have  been  able  to  observe  such  S  curvsHS 

in  plots  of  sales  or  market  share  as  functions  of  the  number  or 

1  3 
the  relative  number  of  outlets 

The  theoretical  arguments  in  favor  of  such  a  S   sha- 
ped relationship  at  the  market  share  level  -  already  introduced  in 
the  form  of  equation  [15)  -  may  also  hold  at  the  transition  proba- 
bility level.  Thus,  at  one  extreme,  if  the  oil  company  has  too  few 
filling-stations  consumers  will  notice  them  too  infrequently  and 
will  often  be  obliged  to  tank  up  at  other  companies  stations; es  a 
consequence,  their  loyalty  will  be  very  low.  As  the  number  of  sta- 
tions increases  consumers  will  be  able  to  tank  up  at  the  company 
petrol  pumps  locatod  in  various  geof.r  aphi  c  a  1  areas  and  their  loyalt 
will  be  enhanced  accordingly.  At  the  other  extreme,  if  the  company 
continues  to  extend  its  distribution  network,  each  new  station  conr 
tructed  wil]  i,::"r  to  attract  consumers  from  the  remaining  hard  cn:- 
of  competitors  customers. 


Estimatinc  thg  pararr.ctors  of  equation  (20)  v/ill  raise  on  p.    isr-ue. 
If  thR  hir.  toricel  data  come  from  a  stable  marknt,  i.e.  our  cutlet 
shares  and  markrt  F^hares'-.  hovv  relatively  little  variatjility.  thei-e 
will  be  severe  nulticnllinee'rity  probl^uos.  As  mnritioned  in  sp':;tion 
3,  the  independent  vr. riahles  in  the  linf;ar  equation  derived  frorr. 

a  first  orrinr  Taylor  t>/:pansiori  are  first  order  derivatives  evalua- 

o 
ted  at  the  current  solution.  Let  the  current  solution  be  a.  =  a.  , 

1     1 


(22) 


3m    [  a  .  =  a. 


^.t-1  °i,t  / 


i.t 


(23) 


^^i.t^^i=^l 


b  .  =  b  . 
1   1 


'i.t-1 


i. 


tween  8rn.   /3a.  and  3m   ./3b, 


An  important  issue  in  nonlinear  estimation  is  fin- 
din^    good  initial  values.  VJe  used  the  follov/inj'  procedure.  In  stead 
state,  the  elasticity  of  market  share  with  respect  to  the  relative 
number  of  outlets  is. 


(24) 


A  .  •  m 
1     c  ,  t  ,  e 


13. 


For    thpsp    rnasons,     a    S       slu^peci    curve    relating,    tran- 
sition   probabilitlDs     to    the    ralative    num'npr    of    retrtil    outlets    '^poms 


(18)  A.     -D.       /(a.     +n.     3 

11  11 


approaches  onr,  according  to  a  S  shaped  pattern 


Similarly  X   is  defined  as 


b  b 

(19)      \       -    D    ^    /     [o       *    D    ^ 

C       C  CO 


with  a   and  b   positive  constant  r,  .  The  s^'^itchinp,  pmb  abilities  are 

b.  b^ 

defined    as    0       =a./(a.     +D.     ),       a.     =a/(a       +D       ). 
ciii  iccc 


Market  share  at  tirr^e  t  is  then 


bi      ^c  ^i  ^c 

^^°^      %.t  -'     ff"i,t  •  °c.t  -  ^  ^^  '     '-\    '    °i.t^-^%  ^  ^c.t 


I.  .  .   +  [a   /  (a   +  D   ^)  1 
i , t- 1       c      c     c. t 


If    D       remains     equal     to     D.        .     and    D       to    0  .     steady 


state    market    share    v/oulri    be, 


(21) 


1,     ^  =a(a.     +D.     ^)/[a.(a       +D       ^}     +     a(a.     +D.     .) 

i.t.e  ci  i,t  ic  c,t  ci  i.t 


To  estimate  initial  values  for  the  parameters,  we  assumed  a  value 


value  for  m 


in  a   =  0.62S  and  b   =  ?.46. 
c  c 

With  these  initio  1  values  TROLL  failcid  to  find  an 
optimal  solution.  Divergence  ocf.  urrori  for  these  and  for  all  other 
initial  values  which  we  tried.  The  SUMT  program  performed  hetter. 
The  main  reason  is  pro  ti  ably  that  the  method  for  minimizinp,  the 
unconstrained  penalty  function  is  the  Nowton-Rephson  method,  a 

second  order  procedure,   whereas  f^larouard  t  '  s  estimation  method 

1  '♦ 
only  uses  first  order  derivatives   .  Furthermore,  all  our  SUMT 

"computations  were  done  in  double  precision.  This  may  be  pai-ticular- 
ly  important  in  view  of  the  factthat  there  is  a  real  multicollinea- 
rity  problem.  The  SUf-IT  estimates  are>  shov;n  in  Table  VI.  The  t  sta- 
tistics ere  very  poor,  as  expected.  Weverthelcss,  the  coefficients 
all  have  the  correct  sign,  and  the  magnitudes  are  reasonable.  After 
all,  the  estimated  values  are  not  too  different  from  our  initial 
subjective  estimates. 


-t  Table  V,  about 


to  an  actual  market  share  of     .  Vc'^le  V/II  shows  the  incremental 
sales  for  adding  a  now  outJet  in  each  of  the  r8i:;lens. 
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sert  Tabl-  VI  I  ab?ut  her 


4.3.  _no_do_l_  £P£lj_f^a.t  i  Oj2 

TflblRS  \1  and  VIT  shov-;  rpspcctively  how  the  e>;ponf?n- 
tinl  and  lop,  istic  morials  can  be  applied  to  compare  in  crp  mental  sa- 
les from  invGstmRnt  in  outlets  in  various  regions. 

The  ultimate^  measure  of  performance  is  the  nompari- 

1  5 
son  of  rofional  profits  rather  tl-ian  sales   .  The  cost  structure 

may  differ  from  cne  rep,  ion  to  another.  For  example,  trans  [lortation 

costs,  and  cost  of  purchasing,  land  will  vary  a  i:  cross  regions. 

The  optimization  model  (R]  proposed  by  H-F  could 


be  adjusted  for  the  exponential  transition  [probability  functions  o- 

section  4.1.  We  replace  q    •  ri"''   in  the  objertive  function  by 

0''^  •  ^i  t  •  """^^  constraint  q '^  ^  =  k^-Q-^^  /  (d^^  +  y  '  ^'[^  ]     is  repla- 

16* 
c  e  d  b  y    : 


i  ,  t 


c   c.t   _j 


c   c.t 


i,t-1 


or  if  cur  interest  is  only  in  steady  state  results  without  conceri 
for  the  transient  behavior. 


-aJ.D^  , 


''i  "i.t 


J.n- 


I'^th  the  loeistic   model,  complications  would  ai-i;:t 
in  the  optimiz'tlon,  because  the  transition  probability  functions 


77 


are  convex  for  some  ronrc  of  D.  [or  D  ),  and  concave  clEewhRrn. 

1        c 


NeverthG]o3S  the  model  will  remain  useful.  Instead 
of  trying'  to  find  the  optima]  allocation  over  time,  the  model  could 
be  applied  to  civaluate  various  outlet  expansion  plans. 

5.  CONCLUSION 


In  this  paper  we  have  presented  a  detailed  study 
of  the  ag^eref, ate  retail  outlet  location  model  deve]opcd  hy  Hartung 
and  Fisher.  V.'e  found  that  their  E^'^eral  procedure  was  sound,  but 
that  the  specifics  suffered  from  a  variety  of  problems  rep, ordinp, 
the  estimation  of  thij  model  parameters,  and  the  i-obustness  of  the 
response  functions.  Rem.edial  action  was  proposed,  and  was  applied 
to  a  brand  of  gasoline  in  a  European  country.        .^^ 
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TABLE    III 


Steady  State  Estimation  of  Exponential  rodel 
for  brand  i 


Coefficient  Value        t  statistic 


-  26.241 


3.25 


25 


.  700 


F(2  >33)  =  3i 


DW  =  0.49 


SUMT 


TABLE    IV 


Exponential  Modnl  for  Brand  i 


26, 


Coef f Icinnt 


Statistic 


F{?  .32) 


59. OG 


2.  75 


.70 


68 


0.  73 


2.76 
7  .  63 


70 


69 


37.3 


TABLE  VI 


SUM    Estimates    for    Logistic    ModGl 


?'3 


CoRfficient  Value 


t  statisti 


8394 


.019 


2.663 


.061 
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FOOTNOTES 


Based  on  discussions  with  various  major  oil  companies.  The  cGnciir;t 
of  linked  hierarchical  models  has  hpen  examined  hy  CrL.vJstcn  and 
Scott-Horton  \  1  ]  ,     and  has  been  applied  in  the  area  of  operatior-ril 
planning  and  control  by  Green  FB],  New  son  [15],  and  Shwimer  T'^Bl. 
A  possible  application  in  the  public  sector  fias  been  proposer'  hy 
Hausman  and  Naert  [6]. 

Our  notation  will  differ  from  H-F. 

See  tl-F  \  7  ,     p.  R  234,  footnote  ?].  The  quote  is  given  in  our  nn- 
tation . 

At  least  within  the  limits  discussed  above. 

TROLL  is  an  interactive  estimation  and  simulation  p-jckage  develop 'ic 
by  the  Computer  Research  Center,  National  Bureau  of  Economic  Re- 
search, Inc.  For  a  description  of  the  estimation  riethod  see  [3]* 

For  a  further  discussicn  on  this  method,  '.-.'o  may  refer  for  cxa'^'.pl" 
to  Draper  and  Smith  f?,  pp.  267  et  sq.l  and  GOIDFELD  and  QUA^DT 
[5.  pp.  49-57] . 


'''    Provided  that  the  linearized  form  (11)  of  the  model  is  valid 
6,  the  final  estimate  of  6. 


^    A  test  of  hypothesis  conducted  hy  LAhLilN  Fill. 

^    For  brand  a,  38  observations  v.' ere  available. 

^    At  this  stof^c  an  additional  qualification  ought  to  be  made  about 
[■l-F  estimation  procedure.  The  r.  imilarity  of  the  results  oritainsd 
from  regression  analysis  apnlied  on  equations  (fi)  and  (5)  as  oppo- 
sed to  the  relative  dissimilarity  observed  betvveon  those  derived 
from  equations  (7)  and  (5)  should  cause  no  surprise,  since  equa- 
tion [ 7 )     contains   m. ,  a  stochastic  variable,  on  both  sides  [as 
evidenced  by  equation  (fi)].  fie  nee  we  sfiould  expect  biased  estimates 
of  a    and  6  from  small  samples  ;  furtherrnorc;  in  this  case  a  and  6 
are  inccnslstent  since  the  residuals  indicate  a  strong  autocorrela- 
tion among  the  error  terms. 

All  this  adds  up  to  pointing  out  the  bias  created 
in  estimatir'g  K.  and  k„  via  equation  (7). 

Notice  that  no  error, term  was  introduced  in  equatior 
(5], had  one  be  added  the  linearization  of  (5)  would  have  been  impos- 
sible. 

^    For  additional  discussion  on  li.  nearizinc,  such  nonlinear  models,  set 
Naert  and  Bultez  [14  1. 

^    For  a  theoretical  exposition  of  SUMT,  see  Fiacco  and  Mc  CormicK 
[41.  The  computer  program  is  described  in  MylancJer  et  al.  [13]. 
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'■  ^  Many  factors  making  the  S  shaped  curve  plausible  in  this  industry 

are  reported  by  Kotler  [9,  pp.  96-971. 

^ '*    See  Wilde  and  Bcightler  fl7,  pp.  22-24]  for  a  discussion  of  the 
Newton-Raphson  method.  The  SUMT  routine  has  various  options  with 
regard  to  the  technique  for  minimizing  the  penalty  function.  One 
of  these  is  the  Newton-Raphson  method. 

^^    Assuming  that  orofit  maximization  is  the  objective. 

^^    Different  types  of  regions,  e.g.  metropolitan, suburban,  rural,  m.igh 
have  dif f ei 
by  regions 
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